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The continuity properties of the radial distribution function g{r) and its close relative the cavity 
function y{r) = e'^'^'^'^^''^'^ g{r) axe studied in the context of the Percus-Yevick (PY) integral equation 
for 3D square-well fluids. The cases corresponding to a well width, (A — l)cr, equal to a fraction 
of the diameter of the hard core, cr/m, with m= 1,2,3 have been considered. In these cases, it is 
proved that the function y{r) and its first derivative are everywhere continuous but eventually the 
derivative of some order becomes discontinuous at the points (n + l)cr/m, n = 0, 1, . . .. The order 
of continuity (the highest order derivative of y(r) being continuous at a given point), k„, is found 
to be Kn ~ n in the first case (m = 1) and ~ 2n in the other two cases (m = 2, 3), for n 3> 1. 
Moreover, derivatives of y{r) up to third order are continuous at r = ct and r = \a for A = 3/2 and 
A = 4/3 but only the first derivative is continuous for A = 2. This can be understood as a non-linear 
resonance effect. 
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I. INTRODUCTION 

Simple models have played an important role in the development of liquid state theory, both as approximations and 
as useful reference systems in perturbation schemes. The simplest model with a repulsive hard core and an attractive 
well is the square-well (SW) fluid. The SW interaction potential is 

oo, r < a 

-e, (J <r < \(7 (1) 
0, r > Act, 

where a is the diameter of the hard core, e is the well depth, and (A — l)cr is the well width. The equilibrium properties 
of the fluid depend on the values of three dimensionless parameters: the fraction of volume occupied by the spheres 
rj — 7r/9(T^/6 (p being the number density), the reduced temperature T* — kBT/e (T being the temperature and ks 
being the Boltzmann constant), and the width parameter A. These properties can be derived from a more fundamental 
quantity, the so-called radial distribution function g{r) The quantity 'i-Kr'^ g{r)dr /V is the probability that the 

centers of two spherical atoms in the liquid are separated by a distance between r and r -I- dr, V being the volume of 
the system. In the absence of particle interactions (ideal gases) there is no structure in the fluid and g{r) = 1. This 
is also true for any realistic potential in the limit r — > c» as we must have (/)(r) — > in that limit. It is also evident 
that the centers of two particles in the SW fluid cannot be nearer than a due to the hard core, so g{r) = for r < cr. 
Many closed integral equations have been proposed for this function (YBG, HNC, Percus-Yevick, . . . ), but the most 
popular is, perhaps, the Percus-Yevick (PY) equation. In the search for these approximations it has been found useful 
to define an auxiliary function, the direct correlation function c(r), through the Ornstein-Zernike (OZ) relationship 

h{\v2 - ril) = c(|r2 -r,\)+pj d^rgcdn - rgDMIra - rsl) (2) 

where h{r) = g{r) — 1. In this equation and the following ones we are assuming that the interaction potential is 
spherically symmetric and the structure functions, consequently, depends only on the distance between the particle 
centers. The physical meaning of the direct correlation function is clear; it accounts only for the correlation effects 
arising from a direct interaction between the particles in the fluid. Nevertheless, the total correlation function, h{r), 
comprises also the correlations propagated by intermediate particles. The direct correlation function is expressed in 
terms of the so-called bridge functionals, which constitute an infinite sum of diagrams Integral equations are 
obtained by closing the OZ relation with an approximation for c(r) based upon a sum of some set of these diagrams. 
By choosing an appropiate set of diagrams, Percus and Yevick showed that the direct correlation function is 
approximately given by 

c(r) = g{r) \l - e^('-)/'=«^l (3) 
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According to (||) the range of c(r) is equal to that of the interaction potential (f>{r) as c{r) is zero in those regions 
where the interaction potential vanishes. The PY equation is given by the substitution of (|^ into the OZ relation (||) 
and it takes the form |^ 

y(|r2 ^ri\) = A + pJ dh^ /(jn - raj) [1 + /(jra - r^l)] y(|ri - rgDydrg - raj) (4) 

where y{r) = e'^'^'^^/^'^'^ g{r) is the cavity function, /(r) = ^-'t>(r)/kBT _ ^ is the Mayer function and A is a constant 
defined as 

/•oo 

A = l- p I drf{r)y{r)4TTr'^ (5) 
Jo 

The PY equation is exactly solvable for the hard-sphere (HS) fluid and the sticky- hard-sphere (SHS) fluid [||. 
Both models can be considered as special cases of the more general SW interaction: the HS potential is recovered if 
we take the limit e ^ (i.e., T* — > oo) or A ^ 1; the SHS fluid is defined by taking the limits A ^ 1 and e ^ oo (i.e., 
T* 0) simultaneously, while keeping the parameter = 12 (l — A~^) e^/"^ constant. 

Nevertheless, neither the PY equation nor any other integral equation for fluids has ever been analytically solved 
for the SW interaction. In 1977, Sharma and Sharma Q proposed a mean spherical approximation which provides 
an analytical expression for the structure factor but it is not consistent with the hard core exclusion constraint. By 
the same time, Nezbeda [|| proposed a polynomial approximate solution for y{r) valid in the limit A — 1 <C 1. This 
solution was based on the continuity of the first and the second derivatives of y{r) at r = a. More recently, Yuste 
and Santos j^] derived a simple approximate expression for the Laplace transform of rg{r) for the SW fluid. This 
derivation is based upon simple physical requirements (flniteness of the radial distribution function at r = (t+ and 
finite isothermal compressibility) and the continuity of y{r) at r = A. 

The aim of this paper is the analysis of the continuity properties of the function y{r) satisfying the PY equation 
for the three-dimensional SW fluid. First, we must identify the points at the borders of the regions where y{r) is 
an analytic function. For A = 2 it is clear that these points are given by r = cr, 2a, . . ., as also happens for the HS 
potential [Q. The radial distribution function is then divided in analytic pieces "00 (^): ■ ■ ■ ^-nd the PY equation 

(^ is written as a system of non-linear integral equations for them. In these equations the Heaviside step function 
that enters in Eq. (m through the Mayer function no longer appears. The continuity properties are derived from this 
system. The cases A = 3/2 and A = 4/3 are more cumbersome as the intervals cited above are divided into two and 
three equal parts, respectively, and the number of equations is obviously larger. More general cases with 1 < A < 4/3 
are increasingly more difficult to manage since the number of intervals becomes larger and larger as the well width, 
A — 1, becomes smaller. We will not deal with those cases in this paper. 

The paper is organized as follows. In Sec. || the PY equation is set in a form that fits better our purpose and the 
continuity of y{r) and its first derivative for r = Act is suggested from the simulation results of Henderson et al. [nOl. 



The appropiate system of integral equations for the case A = 2 is written in Sec. HI; the general derivatives are then 



related with the derivatives of lower order and the continuity conditions are derived recursively. A similar procedure 



is used in Sec. [V for the cases A — 3/2 and A — 4/3. The conditions for the HS fluid are given in an Appendix. The 



paper ends with some remarks on the relevance of these conditions for the proposal of approximations. 

II. DEFINITIONS AND BASIC EQUATIONS 

In the applications of the PY equation to simple potentials is often useful to define a new structure function 
i}{r) — ry{r), which satisfies the same continuity conditions as y{r). After the change of variables t = \r2 — ri|/c, 
cc = |ri — ral/a and y — — \ / a one can easily check that Eq. (^) takes the form j^] 

m=At + 12Tj dxf{x)^{x) dy[l + f{y)]^{y) (6) 

JO J\t~x\ 

where 

POO 

A = 1 - 24ry / dxf{x)x4}{x) (7) 







It is remarkable that this reduction is only possible in three dimensions (excepting the one-dimensional case, where 
no reduction is necessary) and that the same change of variables in two dimensions yields a significatively more 
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complex integral equation due to the lack of rotational simmetry around the vector r2 — ri. The Mayer function 
f{x) = exp [— 0(x)/fcBT] — 1 for the SW interaction (|^) is given by 

l + f{x)^H{x-l)[l + aH{X-x)] (8) 

where H{x) is the Heaviside step function and a = exp (c/kBT) — 1. The first derivative of tp{t) obtained from (^) is 
also given here for future reference 




FIG. 1. Monte Carlo simulation results [|lO|| for the structure functions g{r) (squares), y{r) (circles), and i/'(r) (triangles) on 
a SW fluid with A = 3/2, T* = 2/3 and rj = 27r/15. Note that g{r) and y{r) overlap for r > Act. 



ip' (t) ^ A+ 12?7 / dxf{x)ilj{x) [1 + f{x + t)] iIj{x + t) 
Jo 

-12r/ / dxf{x)^{x)[l + f{t^x)]^{t-x) (9) 
Jo 

/oc 
dxf{x)^{x) [1 + f{x - t)] ^{x - t) 

In Fig. |l]the simulation results of Henderson et al. 0] for g{r), y{r) and il^{r) for a SW fluid with A = 3/2, T* = 2/3 
and rj — 27r/15 are shown. The radial distribution function is not continuous at r = Act but both y{r) and i^{r) are 
clearly continuous functions at that point. Simulation results also suggest the continuity of the first derivative but 
nothing can be said about the second or higher order derivatives as these are very difficult to measure in a simulation, 
where only results for widely spaced values of r are provided. 



III. THE SQUARE- WELL FLUID WITH A = 2 

In this case the well width coincides with the diameter of the hard core and the function that satisfies the 
corresponding PY equation (H) is conveniently defined in a piecewise fashion as follows 



' "00(7-) 

V'i(r) 



< r < 1 

1 < r < 2 



(10) 



7/'„(r) n <r < n+1 
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where we assume, to be confirmed later, that ipnif), n = 0, 1, . . . are analytic functions in their intervals of definition. 
This will be justified by expressing (||) as a system of integral equations for these functions where the Heaviside step 
function in (||) no longer appears. Thus, any derivative of Tpn{r), n = 0, 1, ... is related to derivatives of lower order 
of the same set of functions. It can be expected that at a given point t = n+1, n = 0, 1,... the functions ipnif) 
and ipn+iif) do not match perfectly and this gives rise to a difference between the derivatives of some order of these 
functions evaluated at that point. We define the symbols 

Ai'^ =^it\n+l)~^^l,in+l) n-0,1,..., fc = 0,l,... (11) 

as the jump in the derivative of order k evaluated at i = n + 1. The continuity of ip{t) is a well known fact [^,^, 

so that we have A^"-* = 0, n ~ 0, 1, . . .. The first derivative is given in (||) in terms of integrals over the functions 
■0(0 and /(t), which in turn depends on the Heaviside step function. Starting from (^ and (^ the following set of 
equations for ■0„(t), n = 0, 1, . . . is found 

%{t) ^ A - 12i]{a + 1) [ dxtpo{x)ipi{x + t) 
Ji-t 

2 p2-t 

dx'i(ji{x)'ipi{x — t) + / dx'ipi{x)ipi{x + t) 
+t Ji 



+12r]a{a + l) |^ 
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1277a / dxtpi{x)^p2ix + 1) (12) 



ip[{t) = A - 12r] ^^J dx'ipoix)iJi{x + t) 

+ / dxipo{x)'ip2ix + t) - / dx^po{x)^pl{t — x) 

J2-t Jo 

+ 12770; |y dxipi{x)'i(j2{x + t) + J dxipi{x)^3{x + t) 

2-t (.t-1 

dxipo{x)'ipi{x + t) + / dx'ipo{x)4'i{t — x) > (13) 

Jo } 

■ip'2{t) = A - 12i] [ f dxipo{x)ip2{x + t) + f dxtlJo{x)^p3{x + t) 
Uo Js-t 

1-1-2 r-l 

dxijjo{x)'4'2{t ~ x) — / dxipQ{x)'4'i{t — x) 

lO Jt-2 

+12770; < / dxil}i{x)'ipz{x + t) ^ j dx'i{ji{x)'ip4{x + t) 

t-1 pi 

dx^jJl{x)'l{Jl{t — x) + / dx'ipo{x)ipi{t — x) 

1 Jt-2 

nt-1 

-Urjo^ j dx^i{x)ii^i{t - x) (14) 
tp^{t) ^ A - 12t] I / dx'ipo{x)tpn{x + t) + / dx^po{x)i{jn+i{x + t) 

iJo Jn+l-t 

pt — n pi 

dxipo{x)lpn(t - X) - / dxtpQ{x)'4)n-l{t - x) 
10 Jt-n 

r pn+2-t p2 

-12770 < / dx7/'i(a:)7/'„+i(x + t) + / dxij:i(x)iljn+2{x + t) 

yJl Jn+2-t 
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t-n+1 

dx^l{x)tpn-l{t - X) - I dxipi{x)lpn-2{t - X) 

t-n+1 



— 12r]a / dx^pi{x)^i{t — x)6n,3 , n>3 (15) 

Jt-2 

With the PY equation written in the form of this system it is easy to show that A^^^ = 0, n = 0, 1, ... as a consequence 
of the cancelation of some of the integral terms at the border points. The continuity of tjj{t) and its first derivative is 
also true in the cases A — 3/2, A = 4/3, and in the HS fluid and it is possibly a general property of the PY equation. 
The derivatives of order k, k > 2 of '(pi{t), '02(i), etc. . . are obtained from the system of equations (p^)-(p^), yielding 



+a{a + 1)4'^ (2 - t)7/;f "'"^'^ (2) - a^'^ (2 - t)V;f"'"'^ (2) } + [Int] , (16) 
4''\t) = 12ryX] (-IF {4'^ (2 - t)A^f'^ + (1 + c.)i~iy4'\t - 1)4'"'"^'^ (1) 

-a4'\3 - t)Ai'-'-'^ + ai>l'\2 - t)4'-'-'\2)} + [Int] (17) 
V'f (0 = l2rj''£{i-iy4^\3-t)At'-^^-4^\t^2)At'-'^-a{^^^^ 

3=0 

-ai^l^\t - 2)i^^^-^-'\2) - a{a + l)^p['\t - l)vf ^'^^'^Cl)} + [Int] (18) 
^«(t) = 12rjJ2{i-^yy^i'\n+l~t)Ait'-'->-y,i^\t-n)^^^^^ 

-a(-l)V?' (« + 2 - AiVi'-^' + (t + I - ri)A^:t'^ } 

fe-2 

+12r;a2<5„_3 ^ (t - 2)4'''^'''^ (2) + [Int] , n > 3 (19) 

3=0 

where the terms denoted by [Int] include the sum of several integrals over the functions tpnit), n = 0,1, . . . and their 
derivatives. These terms are always continuous and their explicit expressions are not required for the calculation of 
the derivative jumps, so they will not be quoted here. If we take into account the continuity of '4'{t) and its first 
derivative, as well as the condition '0(0) — 0; the first nonzero derivative jump at the border points t = 1,2, .. . is 
readily derived from Eqs. (16)-(|l9|): 

A[,^^ = A^^^ = -24rya(a + 1)V'(1)V(2) (20) 

Ap' = 12r;(a + l)MV'(l)]' (21) 

Af ) = 12r,a^ mf (22) 

Agr^ - -12ry(« + l)^(l)Agr-f ^ + l2mm^t-2^ , n = 2, 3, . . . (23) 

A^^r^i = 12,7a0(2)A(^f ) , n = 2,3,... (24) 

Therefore, the order of continuity (i.e., the highest order derivative of being continuous) is k„ = 1 for n < 3, 
Kn — n — 1 for n = 4, 6, . . ., and k„ = n — 2 for n = 5,7,.... The presence of discontinuities in the potential at 
r — (J and r — 2(j gives rise to discontinuities in the second derivative of the structure function 'ijj{r) not only at those 
points but also at their neighbors r = 3a, 4a and, surprisingly, these discontinuities are of the same order as those at 
the hard core and square- well borders. At larger distances the structure function becomes more and more continuous 
and the order of continuity grows linearly. The derivative jumps in (p3|) and (p3) may, indeed, increase exponentially 
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at certain physical conditions but this is balanced with the factorial in a Taylor expansion. These results are then 
compatible with the asymptotic limit g{r) = y{r) 1, or equivalently ijjir) —^r^asr—^ oo. It can also be noticed 
that the first nonzero derivative jump in Eqs. (|20|)-(p^) depends only on the packing fraction i], the temperature 
(through Of), and the values of the function ip{r) at the potential discontinuity points. The order of continuity of the 
function ipir) and the cavity function coincide by definition. The derivative jumps of the second derivative of y{r) at 
r = (7 and r = 2(7 are given by 



r=2(T+ 



dr^ 



48?7 



'.{a + l)y{a)y{2a) 



= -%a + ir[y{af 



(25) 
(26) 



r=2cr- 



which are a consequence of Eqs. (|^) and 1x2 



IV. THE SQUARE- WELL FLUIDS WITH A = 3/2 AND A = 4/3 



In the case A = 3/2 the function ip{t) is also piecewise but we must distinguish between the left and the right halves 
of the intervals n < t < n + 1, n — 0,1, . . .: 



< t < i 
h <t <1 



r Mt) 

■ipnii) n+^<t<n+l 



n < t < n 



(27) 



Consequently, the derivative jumps at the points t = n + 1/2 and t = n + 1, n = 0, 1, 2, . . . are denoted by two sets of 
symbols, A^'^'' and A^i'^ , respectively. These symbols are defined as follows 



A(f)=?f (n+l)-^,i^:^,(n+l) 



(28) 
(29) 



As in the previous case the functions '0(t) and its first derivative are everywhere continuous and these symbols take 

the value zero for fc = and fc = 1, n = 0, 1, . . .. The equations for ip^n\t) and V'n are obtained from (^) and (H) 
after the elimination of the Heaviside functions. A straightforward but lengthy calculation that runs in parallel to 
that of Sec. Ill leads to the following results for the first nonzero symbols at every border point 



A 



(4 



Ai- 



(2 



A 



(2 



^(2n-2 



A 



A 



(2 



(4 



^(2n-2 



(12r;)2(l + a)^(l) {(1 + af m)f + [^(3/2)]'} 

12r,{l + af[M)f 
12r;«2 [^-(3/2)]" 

-12^ {^(l)Ail7^) - aV;(3/2)Ail7^)} , n > 3 

A^^^ = -2477a(l + a)V'(l)V'(3/2) 

-(12,7)2a^(3/2) [a^ [^(3/2)]^ + 2(1 + af [M)f] 

-12,? {(1 + a)V'(l)Ai^T" - «V'(3/2)AiV)} , n > 3 



(30) 

(31) 
(32) 
(33) 

(34) 
(35) 

(36) 



In this case we will denote the order of continuity at t — n + 1 by n = 0, 1, . . . and the corresponding one at 
t = n + 1/2 by k„, n = 0, 1, . . .. According to (|30|)-(^) we conclude that kq = ki = 3, ki = ^2 = kq = '^2 = 1, 
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and K„ = K„ = 2n — 3 for n > 3. It is a remarkable fact that the cavity function y{r) is continuous up to the 
third order derivative at those points where the interaction potential exhibits discontinuities, r = a and r = 3a/2. 
Nevertheless, only the first derivative is continuous at the neighbour points r = cr/2, r — 2cr, r = 5a/2 and r = 3a. 
The situation is rather different in the case studied in the previous Section. There, the discontinuity on the square- well 
is located precisely at a point where the hard core induces a discontinuity jump on the structure function. This special 
disposition of the potential jumps reinforces the discontinuity of the structure functions, y{r) and tp{r), at the lattice 
sites r /a — 1,2, ... on a kind of "resonance" effect. This is also observed at large distances from the hard core as the 
highest order derivative exhibiting no jumps grows as r in the SW fluid with A = 2 and as 2r in the SW fluid with 
A = 3/2, where the resonance condition is broken. 

In the case A = 4/3 we can find discontinuities at the points t = n + 1/3, t = n + 2/3, and t = n + 1, n = 0, 1, 2, . . .. 
In a similar way to the definitions in (O), (psl), and (poh the symbols A^^^^ A^*^^ and A^''^ are introduced as follows 



Ai'^) 



Al'=) 



Al'=) 




t=n+l/3-e 



t=n+2/3-e 



t=n+l- 



d'^ipit) 



t=n+l/3+e 



t=n+2/3+e 



dtk 



t=n+l+e 



(37) 



(38) 



(39) 



where n = 0, 1, . . ., fc = 0, 1, . . .. The calculation of these symbols is completely analogous to that of the previous 
cases save for the existence of three independent analytic functions on the intervals {n,n + 1), n = 0, 1, . . . and the 
details will not be given. The results for the nonzero symbols with the lowest value of k for n = 0, 1, 2 are 



A 



A 



A 



A 



A 



A 



A 



A 



(2 



(4 



(4 



(4 



(4 



(2 



(2 



(4 



A 



(2) 



-24rya(l + a)?A(l)V'(4/3) 



~{12rj)'a{l + af [^.(1)]' ^(4/3) 

(12,y)2(l + a)V^(l) {(1 + m)f + 2a' [^(4/3)]^} 

-(12,7)2aV^(4/3) {a' [^(4/3)]' + 2(1 + [m?} 
(12r7)W(l + a)^(l)[V'(4/3)]' 



12?7(1 



[V^i)]' 



12rja' [^-(4/3)]^ 
-127^(1 + a)Af^ 



(40) 
(41) 
(42) 

(43) 

(44) 
(45) 
(46) 
(47) 



Between these conditions and (K(t)-(B6[) corresponding to the case A — 3/2 we notice some similarities. In particular, 

13 5 

the order of continuity at the points t ~ — , 1, — ,2, — and 3 for A — 3/2 is the same that the order of continuity 
14 7 8 

at the points t = I7 3'^' 3 ^"^^ 3 '^'^^ ^ ~ ^13. Moreover, the first nonzero derivative jump has the same form 
in both cases if we write them in terms of i)(X) and '^(X) (See Eqs. ( ^4| ) and (40), ( |30| ) and (^), etc. . . ). These 
considerations suffice us to conjecture that the order of continuity at the points t = A — 1,2,A-|- 1,2A is fc = 1 and 
the order of continuity at t = 1, A is A: = 3 for any value of A in the interval 1 < A < 2. The special "resonant" case 
A = 2 is excluded from this rule. The first nonzero derivative jumps at these points in terms of the cavity function 
are expected to be given by 



dr'^ 



(A±l)o-+ 

dr'^ 



dr^ 



dr^ 



AtT+ 



(A±l)(j- 

d^ 



(A ± 1)0-2 "^ + a)y{a)y{\a) 



(12??)= 



(ISry)^ 



(1 + a)y{a)[[{\ + a)y{a)f + 2 [aXyiXaf] 



ay{\a) { [a\y{\o)f + 2 [(1 + a)y{o) f] 



(48) 

(49) 
(50) 
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2<j+ 



2A(T+ 



_ fy 
dr^ 

dr"^ 



i\o- 



-§[(l + a)2/(a)]' 



{oiy{\a)f 



(51) 
(52) 



which are a consequence of the generahzation of (|30|)~(|32|), (|3j) and ( |35| ) or the corresponding conditions for A = 4/3. 
Results for the order of continuity versus r j a for the SW fluids we have analyzed and the hard-sphere case discussed 
in the Appendix are summarized in Fig. ^. The main feature of these series of values is their linear (HS) or staircase 
(SW with A = 2, 3/2 and 4/3) behaviour beyond a given distance which depends upon the well width. This is a simple 
mathematical consequence of the explicit convolution nature of the PY equation. It can be shown that a solution 
of this kind of equations with order of continuity p at r = a and p at r = 6 is necessarily a function with order of 
continuity p + p +3 at r = a + 6. This rule is sufficient to explain the trend for r > 3cr on every case plotted in Fig. ||. 
Discontinuities of the structure functions derivatives appear at r — a and r = Act as a consequence of the singularities 
of the interaction potential, (p^r) , and are then propagated in a way we are not familiar with in the domain of linear 
physics. The cases analyzed in detail in this paper give us also some insight on the interplay of the relative positions 
of the discontinuites of the potential and we could suggest an explanation for the qualitatively different behaviours 
of the cases 1 < A < 2 and A = 2 showed in Fig. |2| The order of continuity appears to be generally lower in the 
latter case and we can attribute that to an interference of the singularities of the Mayer function because for A = 2 
the square-well border is placed precisely on a point 
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FIG. 2. Order of continuity of y{r) versus r/a for the SW fluid with A = 2 (squares), A = 3/2 (triangles), and A = 4/3 
(dots), and the HS fluid (circles). The order of continuity is assumed to be infinity if not plotted. 



where the singularities induced by the hard core are propagated to. Conjectures about what will occur with arbitrary 
values of A are discussed in the next section. 



V. CONCLUDING REMARKS 



In this paper the analytic properties of the structure function of square-well (SW) fluids have been studied. The 
SW potential is a simple but useful model of more realistic interactions as it includes not only the volume exclusion 
effect of the hard core but also a finite ran ge a ttractive part. It has been used as a model for fluids, systems of colloidal 
particles, microemulsions and micelles 1 11 - 15|] . Throughout this paper, it has been assumed that the structure function 
ijj(r) = ry{r) satisfies the well known nonlinear Percus-Yevick (PY) integral equation. Despite the apparent simplicity 
of the SW potential, no exact solution of this equation has ever been found. The reduction of the PY equation to a 
system of nonlinear integral equations unveils the degree of difficulty of the problem. The structure function, 'ip{rY 
must be defined by analytical pieces and this is unsuitable for a purely local treatment, except in the hard-sphere 0] 
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and sticky-hard-sphere Q models, where an expression for the Laplace transform of rg{r) was found. 
Although we have restricted ourselves to the cases A = 2, A = 3/2 and A = 4/3, the general pattern for any value A 
seems clear. If A is an integer we have discontinuities only at r/a — 1,2,.. ., the "resonance" condition is fulfilled and 
the order of continuity grows as r for large r. If A is not an integer we must expect the appearance of discontinuities 
at those points which can be expressed as a linear combination of the fundamental scale lengths a (the hard core 
diameter) and (A — l)a (the SW width), na + m{X — l)t7, where n and m are integers. If A is a rational number there 
is some point that can be obtained by using two different pairs of integers (n, m) and the location of the discontinuity 
points repeats periodically after that. That is not the case for irrational values of A. A direct consequence of the 
choice of a width well (A — l)a incommensurate with the hard core diameter a is that the order of continuity at 
r/a = 1,2, .. . will possibly coincide with the corresponding to the HS fluid because the two discontinuities of the 
interaction potential do not interfere in this case. In the cases we have discussed (A = 3/2 and A = 4/3) there is still 
some kind of "resonance" and this gives rise to an order of continuity lower in two units to that of the HS fluid at 
r/a = 4, 5 . . . (See Fig. ^). However, the conjecture expressed in Eqs. (]4^)-(p2|) for 1 < A < 2 is sufficient for the 
purpose of the search of approximate structure functions as only this interval is considered physically meaningful in 
the applications of the model to simple fluids or colloids. We must also emphasize that these results have been derived 
in the context of the PY equation but, referring only to the key features of the structure functions, it is possible that 
they are more general. An analysis of other integral equations for fluids and the Mayer diagrams not taking into 
account by them would be necessary in order to clarify this point. 

This work was already started in the early days of the theory of liquids. Percus showed Q, in the context of the 
PY equation for hard spheres, that analytic breakings on the structure functions should appear whenever distances 
between particle centers along a chain are modified so to accommodate a new particle. This result explains the 
discontinuites on the cavity function derivatives at r = ct, 2a, 3a, . . .that are listed in the Equations (|A4[)-( A5). A 



combination of formal and heuristic geometric arguments were used later by Stillinger |16[ in order to identify the 
diagrams which cause discontinuities in the derivatives of g{r) for hard spheres up to fourth order in the density. 
The most simple diagrams are the chains but we have also double chains and triply-connected diagrams that are not 
included in the PY approximation. The topological change associated with the separation or the approaching of a 
given pair of spheres in the cluster corresponding to the diagram within a distance which avoids the contact of all the 
particles was found to be the origin of a singularity of the structure function. In that way, the lowest-order double chain 
gives rise to a singularity at r = -\/3cr and the lowest-order triply-connected cluster is responsible for a discontinuity 
at r = y/8/3a Careful Monte Carlo simulations carried out by Seaton and Glandt |18| for the fluid of adhesive 
spheres have confirmed the existence of discontinuities of the radial distribution function at r — y^S/Scr, \/3a and 
2a. The latter is the most striking one and this is taking into account by the PY solution because it is originated 
by the chain diagram with two bonds. The more complex spatial diagrams are responsible for the discontinuities at 
r — ^/8/3a and -\/3(t which are not included in the PY solution for the sticky hard sphere potential Nevertheless, 
these two are less prominent features of g{r) than the discontinuity at r = 2a and the PY solution gives a good overall 
agreement with Monte Carlo simulation results . Stillinger jl6| has even suggested that the full set of singularities 
of g{r) for the hard sphere or hard disk case is dense throughout the entire range < r < oo as a consequence of the 
formation or breaking of bonds on random packings when r (the distance between two given particles in the cluster) 
is slightly varied. 

Analytical approximations proposed by different routes are equally interesting alternatives to the solution of the 
integral equations. These approximations are usually based on the imposition of continuity conditions at the hard 
core and the SW border. For example, the approximation of Nezbeda |Q for a very thin SW interaction is based 
upon the continuity of the first and the second derivatives of y{r) a.i r — a. On the other hand, Yuste and Santos 
proposed an approximation on the basis of some physical conditions, among which the continuity of y{r) at r — Act. 
The latter authors are forced to fix one of the parameters of their model at its low density value in order to close 
the system of nonlinear equations that those parameters satisfy. The results derived in Sees. Ill and fV allow the 
imposition of continuity conditions on the first derivative of y{r) at r — \ and this would close the system of equations 
of Yuste and Santos's model in a more convincing way. Similarly, Nezbeda approximation could possibly be improved 
by imposing the continuity of the third derivative at r = ct. Work along this line is currently in progress and will be 
published elsewhere. 



ACKNOWLEDGMENTS 



The author gratefully acknowledges Andres Santos for useful discussions and a critical reading of the manuscript. 
Partial support from the DGICYT (Spain) through Grant no. PB97-1501 and from the Junta de Extremadura (Fondo 
Social Europeo) through Grant no. IPR98C019 are also acknowledged. 



9 



APPENDIX A: THE HARD-SPHERE FLUID 



In this case the structure function is defined in the piecewise form given in ( p^ and the set of equations 
satisfied by the analytic functions ipnit), n — 0,1, . . . is obtained from (|l^)-([l5|) if the hmit a — > is taken (which is 
equivalent to take e ^ or T ^ oo). The derivatives of order k, k >2 has already been given in Eqs. (p^)-(p^. By 
setting a = in these equations we get 



fc-2 



^i'\t) - -12ry^ {-iy4'\^ - + pnt] (Al) 

^f\t) = U^Y^ {{-iy^\^\2 t)A['-'-'^ + 4^")(t - 1)4^-^--'") (1)} + [Int] (A2) 

fc-2 

V'«(i) = 12,7E{(-l)>o'^(" + l-^)^"'"'"'^-^o'^ ' (A3) 

where the terms [Int] again denote the sum of several integrals over the functions ^poit), 'ipi{t), . . . , whose explicit 
expressions are not required as these terms are always continuous at the points of interest t — 1,2,.... The derivative 
jumps aI^\ n = 0, 1, . . ., A; = 0, 1, . . . have been already defined in ( pT|) and the nonzero symbols A^*''' corresponding 
to the lowest value of k for every n can be derived from ( [Al[ )-( |A3| ). The results are 

A^^^ = -1277V^i(l)Af^ - 247^7/^1 (1)4'>(0) = (127^)2 [V-(l)]3 (A4) 

Ai^") = (12ry)" [-^(1)]"+' , n>l (A5) 

where = (1 + ri/2)/{l — 77)^ is the PY exact result [^ for the structure function at the hard core contact point. 
Note that the HS fluid can also be seen as a SW fluid in the limit A ^ 1. 
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